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Abstract 

Reflection asymmetric covariant density functional theory (CDFT) based on the point-coupling 
interaction is established on a two-center harmonic-oscillator basis and applied to investigate the 
Ra puzzle, i.e., the anomalous enhancement of the residual proton-neutron interactions 5V pn for 
Ra isotopes around N = 135. The octupole deformation and shape evolution in the Ra and 
Rn isotopes are examined in the potential energy surfaces in {^2,h) plane by the constrained 
reflection asymmetric calculations. The SV pn values extracted from the double difference of the 
binding energies for Ra isotopes are compared with the data as well as the axial and the triaxial 
calculations. It is found that the octupole deformation is responsible for the Ra puzzle in the 
microscopic CDFT. 
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I. INTRODUCTION 



During the past several decades, the importance of the proton-neutron interaction in 
nuclear structure has been widely recognized It affects many aspects of nuclear struc- 
ture, such as the single-particle energy levels, the nuclear shape transition, the onset of 
deformation, the shell closure, etc. [l-4]. 

Since nuclear masses embody the interactions of all the nucleons, it is possible to isolate 
and extract the interaction of specific kind of nucleons from the masses. In particular, the 
average interaction strength 5V pn between the last protons and the last neutrons in even- 
even and even-Z, odd-iV nuclei can be extracted by the double difference of the binding 
energies {5, 6] as 



5V p e :(Z, N) =-{[B(Z, N) - B(Z, N-2)}- 

[B(Z-2,N)-B(Z-2,N-2)}}, 

5V™(Z, N) = 1 -{[B{Z, N) - B(Z, N - 1)]- 

[B(Z - 2, N) - B(Z — 2,N— 1)]}. 



(la) 



(lb) 



With the atomic mass evaluation published in 2003 (AME03) a systematic investigation 



of 5V pn values throughout the mass surface was performed in Refs. [8M1 1| . It was found that 
the results in regions of strong shell closures and in regions where shape transitions occur are 
especially interesting and are able to reflect structural features. Meanwhile, there are also 



many experimental 12j-[l8[ and theoretical 



19 



2l| efforts which were devoted to investigate 



the SV r 



pn- 



In particular, it was noted in Ref. js] that there are anomalous enhancements of 5V pn for 
221 Ra and 223 Ra with TV = 133, 135 deviating from the general trend of 5V pn values, i.e., the 



so-called "Ra puzzle" js]. Later on, a precise Penning-trap mass measurement on 



223-229 



Rn 



has provided clear evidence of the existence of Ra puzzle and found that 5V pn of the odd-iV 



Ra isotopes shows a well developed peak around iV = 135 which terminates at N = 139 [16]. 

It is speculated that the Ra puzzle is associated with the softness of well-known octupole 
deformation in this region jol, [la] . The octupole correlation is due to the interaction between 
orbital pairs with Al = 3 and Aj = 3 around the Fermi surface. For the nuclei around Z = 88 
and N = 134, there exist octupole pairs (w2f 7 / 2 , vrli 13 / 2 ) for protons and (^2g 9 / 2 , ^lj'15/2) 
for neutrons. Therefore, to understand the phenomenon of Ra puzzle, it is necessary to have 



a reliable theory including the reflection asymmetric degree of freedom. 

The covariant density functional theory (CDFT) has achieved great successes in de- 



scribing nuclear properties of both stable and exotic nuclei 
achievements in nuclear magnetic moments 
excitations 
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22H24j . including the recent 
25r|27(, pseudospin symmetry |28l-l33|. low-l ying 



361 ]. magnetic and antimagnetic rotations |37H4l|. collective vibrations 
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461 ]. and so on. Therefore, the CDFT with reflection asymmetry is an appropriate choice 



to investigate the Ra puzzle. In most of the successful versions of CDFT in nuclei, the 



Fock terms 



471 ] are not included explicitly, which leads to the relativistic mean-field (RMF) 



theory and forms the basis of its widespread applicability at present. 

The reflection asymmetric relativistic mean-field (RAS-RMF) theory with meson-exchange 



interaction has been independently developed in Ref. 48[ on grid and in Ref. 49] on a two- 
center harmonic-oscillator (TCHO) basis 50j. This model has been successfully applied in 
the description of the ground-state properties of 226 Ra 4^] and the shape evolution of Sm 51] 



and Th [52] isotopes. In Ref. 53|, |54j, a RAS-RMF theory using both the meson-exchange 



and the point-coupling interactions with the triaxial degree of freedom was developed on 
the conventional harmonic-oscillator basis. 

Recently, CDFT with the point-coupling interaction has attracted more and more atten- 
tions due to its simple applicability in being extended beyond the mean-field approxima- 



tion 



55H57J. In this paper, the reflection asymmetric RMF theory with the point-coupling 



interaction (RAS-RMF-PC) is developed on a TCHO basis. With the RAS-RMF-PC thus 
implanted, the potential energy surfaces in (fa, $3) plane for Ra and Rn isotopes will be 
studied and the relationship between Ra puzzle and octupole deformation will be examined. 



II. THEORETICAL FRAMEWORK 

The starting point of the RAS-RMF-PC model is an effective Lagrangian density with 
the zero-range point-coupling interaction between nucleons. By means of the conventional 
variation principle, the Dirac equation for nucleons can be obtained 

[-ia ■ V + P lt y» + p(M + S)]Mr) = e k Mr), (2) 

where 

S(r) = a s ps + PsP s + lsp s + S s Ap s , (3) 
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+T 3 a T vjTv + T 3 6 TV Aj£ v + eA». (4) 

More details can be found in Refs. 55l-l57l|. The basis expansion method is widely used to 
solve the Dirac equation. For nuclei with reflection asymmetry, the Dirac spinors could be 
expanded in terms of the eigenfunctions of the TCHO potential 

V(r ± ,z) = -Mulrl+{ ' (5) 

2 [ \Mu 2 2 {z - z 2 )\ z > 0, 

where z\ and z 2 (real, positive) are the distances between the centers of the spheroids and 
the intersection plane, and uj\ (u 2 ) are the corresponding oscillator frequencies for z < 
(z > 0) [49j | . The TCHO basis can be completely specified by three parameters: 5 2 , 8 3 and 



49|. 



Az, and their detailed definitions can be found in Ref. 

The binding energy with a given deformation can be obtained by the deformation con- 
strained calculation, i.e., by minimizing 

(H') = (H)+ 1 -C((Q 2 )- f , 2 f, (6) 

where C is a spring constant, \x 2 is the given quadrupole moment, and (Q 2 ) is the expectation 
value of qudrupole moment operator Q 2 = 2r 2 P 2 (cos8). The octupole moment constraint 
can also be applied similarly with Q 3 = 2r 3 P 3 (cos#). By constraining the quadrupole 
moment and octupole moment simultaneously, the total energies in (^ 2 , f3 3 ) plane can be 
obtained. 

III. NUMERICAL DETAILS 

In the present work, all the RMF calculations are performed with the newly proposed 
effective interaction PC-PK1 5^]. The pairing correlations are neglected at this moment. For 
the axial and triaxial calculations, the Dirac equations are solved on the three-dimensional 
isotropic harmonic-oscillator basis. For the RAS-RMF-PC calculations, the Dirac equations 
are solved on the TCHO basis with 5 2 = 0.0, 8 3 = 0.99 and Az 0. By increasing the 
major shell number of the harmonic-oscillator basis from Nf = 16 to 18, the binding energy 
of 212 Ra changes less than 0.01% for axial and triaxial calculations, and less than 0.02% for 
RAS-RMF-PC calculations. Therefore, the major shell number Nf = 16 is adopted in the 
following calculations. 
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FIG. 1. (Color online) The potential energy surfaces for Ra isotopes in ((32,03) plane in the 
reflection asymmetric covariant density functional calculations with PC-PK1. The energy difference 
between neighboring contour lines is 0.5 MeV. The global minima are denoted by solid stars. 

IV. RESULTS AND DISCUSSION 



In order to investigate the shape evolution in the Ra and Rn isotopes with the octupole 
degree of freedom, the potential energy surfaces in (0 2 , 03) plane for 210_229 Ra and 208_227 Rn 
have been calculated in the constrained RAS-RMF-PC theory. As examples, the contour 
plots for 212 Ra, 215 Ra, 218 Ra, 221 Ra, 224 Ra, and 227 Ra are shown in Fig. [Q and the global 
minima therein are denoted by stars. It is shown that the ground states of 212 Ra, 215 Ra, 
and 218 Ra are near spherical without octupole deformation. For 218 Ra, apart from the 
global minimum with 0% = 0, there exists a local minimum with the octupole deformation 
3 = 0.13. The octupole deformation appears in the ground states of 221 Ra, 224 Ra, and 
227 Ra, and increases with the neutron number. 

By analyzing the potential energy surfaces for all the Ra isotopes, one could further see 
that the global minima for 210_219 Ra are slightly prolate with 2 < 0.1. Noted that for 217 Ra, 
218 Ra, and 219 Ra, there appears one local minimum with obvious octupole deformation. For 
220 " 229 Ra, the octupole deformed minimum becomes the global minimum, which means that 
the ground states of these isotopes are octupole deformed. Quantitatively, for 210 ~ 219 Ra, the 
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FIG. 2. (Color online) The residual proton-neutron interactions 5V pn data (solid squares) [la] for 



Q 



odd-iV Ra isotopes in comparison with the calculated values 
(c) CDFT with PC-PK1, as well as HFB-17 mass table (d) 



5y 



Si- 



axial (a), triaxial (b), and octupole 



ground-state deformation parameters are of 0.04 < 02 < 0.1 and /3s = 0. In particular, the 
ground state of the magic nucleus 214 Ra is near spherical (/^ = 0.04). For the ground states 
of 220_229 Ra, both the quadrupole and octupole deformation parameters become remarkable 
with 0.15 < 02 < 0.23 and 0.1 < 03 < 0.18, and increase with the neutron number. 

The shape evolution of Rn isotopes is similar to that of Ra isotopes. The ground states 
of 208_217 Rn have only quadrupole deformation and the /?2 varies in the range from -0.07 
to 0.11. The octupole deformation appears in the ground states of 218-225 Rn, and the /?3 
increases from 0.02 for 218 Rn to 0.08 for 223 Rn but finally drop to 0.02 for 225 Rn. The 
octupole deformation further vanishes in the ground states of 226 Rn and 227 Rn. Compared 
to the Ra isotopes, the Rn isotopes with the same neutron number have smaller octupole 
deformations. 

From the binding energies of neighboring nuclei, the residual proton-neutron interactions 
5V pn for Ra isotopes can be extracted according to Eq.([T]). In Fig. [21 the 5V pp . values extracted 
from RMF models are compared with the empirical values from the data 16( for odd-iV Ra 
isotopes. Here, the results in the axial and the triaxial as well as the reflection asymmetric 
RMF calculations are respectively denoted by open circles, triangles, and solid circles. 

In Fig. [2] (a), the axial RMF calculations reproduce the data well except the data for 
221 Ra, 223 Ra, and 225 Ra, and thus fails in reproducing the peak around N = 135. The 
same conclusion remains even after the triaxiality is considered, as shown in Fig. [2](b) . This 
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FIG. 3. (Color online) Same as Fig. [2l but for the even-even Ra isotopes. 



indicates that the quadrupole deformation and triaxiality are not the reasons for the Ra 
puzzle. 

After including the octupole degree of freedom, as shown in Fig. [2] (c), the peak around 
N = 135 for the 5V pn value is well reproduced. This clearly indicates that the Ra puzzle 
can be well understood with the octupole deformation. It should be pointed out that the 
discrepancies appear at 227 Ra and 229 Ra, which might be attributed to the pairing correlation 
neglected in present RAS-RMF calculations. 



For comparison, the results from the nuclear mass tables HFB-17 58[ are presented in 
Fig. [2] (d) and compared with the data. Again, the HFB-17 results fail to reproduce the 
peak around iV = 135 due to the absence of the octupole degree of freedom. 

Similar to Fig. [2j the residual proton-neutron interactions 5V pn data for even-even Ra 
isotopes are shown in Fig. |3] in comparison with the corresponding calculations. The data 
for light Ra isotopes are well reproduced in both the axial and triaxial calculations. However, 
the discrepancies for SV pn appear for the Ra isotopes from N = 132 to 136. After including 
the octupole deformation, the data are well reproduced expect for 224 Ra and 228 Ra. 

Comparing with Fig. [2]^c), the agreement with the data for the even-even Ra isotopes is 
less impressive than that for the odd-iV isotopes. Further improvement may be achieved 
by taking into account the pairing correlation. In fact, as shown in Fig. [3^d), the results 
from HFB-17 mass table with the pairing correlation could reproduce the data quite well. 
Therefore, in the relativistic framework, both the octupole deformation and the pairing 
correlation might be important to describe the 5V pn values for the even-even Ra isotopes. 
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FIG. 4. (Color online) Neutron single-particle levels for 219_225 Ra in the RS-RMF-PC (left) and 
RAS-RMF-PC (right) calculations. The solid circles denote the levels occupied by the last neutron. 
The levels are labeled by Nilsson-like notations £l[Nn z mi] of the first component. On the top margin 
are the deformations. 

For odd-iV isotopes, the SV pn in Eq. [TJb) could be rewritten as 



with the single-neutron separation energy S n (Z,N) = B(Z,N) — B(Z,N — 1). For odd- 
N Ra isotopes, they are the single-neutron separation energy differences between the Ra 
isotopes and the Rn isotopes. In the RAS-RMF calculations, for Rn isotopes, as their 
octupole deformations are quite small, the S n calculated with and without octupole degree 
of freedom is similar. For Ra isotopes, as their octupole deformations are remarkable, the S n 
calculated with and without octupole degree of freedom is considerably different. Therefore, 
the difference of 5V pn values with or without octupole deformation is determined by the 
corresponding difference of the single-neutron separation energy for Ra isotopes. 

As the single-neutron separation energy is nothing but the Fermi surface in the single- 
particle spectrum with pairing neglected, it is interesting to examine the effect of the oc- 
tupole deformation on the neutron single-particle levels, as shown in Fig. HJ for 219 ~ 225 Ra 
obtained by the RMF calculations with and without reflection symmetry. For 219 Ra, the 
same single-particle level structure is obtained as (3$ = 0. From 220 Ra to 225 Ra, the level 
structure is quite different due to the octupole deformation. In the axial calculations, the /3 2 
deformation increases with the neutron number and the shell gap at magic number A" = 126 




(7) 



S 



in 219 Ra gradually disappears. In RAS-RMF calculations, however, due to the performance 
of octupole deformation the shell gap at N = 126 disappears and a large energy gap at 
N = 132 appears for 220 ~ 225 Ra. 

By switching on the octupole deformation, the Fermi energies 

of22 0-225 Ra 

are bound more 

deeply, which results in larger single-neutron separation energies and provides reasonable 
explanation for the Ra puzzle, as shown in Fig. [2(c). 

V. CONCLUSION 

In summary, the reflection asymmetric CDFT based on the point-coupling interaction 
is established on a two-center harmonic-oscillator basis. The potential energy surfaces in 
(02,03) plane are calculated by the constrained reflection asymmetric calculations to inves- 
tigate the shape evolution in the Ra and Rn isotopes. It is found that the ground states 
are near-spherical for 210_219 Ra and have remarkable octupole deformation for 220_229 Ra. In 
comparison, the Rn isotopes have smaller octupole deformations but similar shape evolutive 
behavior as the Ra isotopes. 

The residual proton-neutron interactions SV pn for Ra isotopes are extracted from the 
double difference of the binding energies of Ra and Rn isotopes and are compared with 
the experimental values as well as the axial and triaxial RMF calculations. It is found 
that the octupole deformation provides a reasonable explanation for the Ra puzzle, i.e., the 
anomalous enhancement of 5V pn for Ra isotopes around iV = 135. 

This explanation for the Ra puzzle by the octupole deformation can be traced back to 
the single-neutron separation energy and the single-particle energy spectrum. The octupole 
deformation will drive the Fermi surface to be bound more deeply for 220_225 Ra, which 
results in larger single-neutron separation energy as well as the appearance of an energy gap 
at N = 132 for 220 ~ 225 Ra. 
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